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I. INTRODUCTION 



Gravitational radiation by non-relativistic and quasi-relativistic systems is low-multipole and can be easily 
calculated using the quadrupole formula of General Relativity with higher multipole corrections. With increasing 
velocities, the contribution from higher multipoles becomes dominant, so one needs another technique. The most 
adequate approach is the method of post-linear expansions which was discussed in early 60-ies most notably 
by Bertotti Q, Bertotti and Plebanski Hi, Havas and Goldberg [1,0 (see also (HHI). We have developed a 
momentum space version of this approach [15j which is applied here to gravitational bremsstrahlung. Although 
technically different, our calculations essentially overlap and agree with those by Thorne and collaborators (la - 
Il8j |. The results also agree with an alternative calculations by Peters (UHH based on the linear perturbation 
theory on Schwarzschild background. They disagree, however, with calculations based on "equivalent gravitons" 
method [22j], and we explain the origin of this disagreement. 

Other approaches to relativistic bremsstrahlung problem are worth to be mentioned. One, suggested by 
D'Eath [23( 1 . is based on replacing the boosted Schwarzschild metric by the impulsive gravitational wave. An- 
other, due to Smarr 24, 25], appeals to calculation of the radiation amplitude in the low-frequency region. This 
overlaps with quantum calculation of the cross-section in the Born approximation [26| . 

To calculate the leading order gravitational radiation in relativistic collisions of particles interacting predomi- 
nantly through the non-gravitational forces it is enough to use the linearized gravity on Minkowski background. 
In the case of gravitational interaction we need at least the next post-linear order. If one interprets the second 
order gravitational potentials in terms of Minkowski space coordinates, one finds that the source of gravita- 
tional radiation becomes non-local due to contribution of the first order gravitational stresses (similarly for 
non-gravitational radiation from particles interacting by gravity). This non-locality leads to destructive inter- 
ference of high frequency part of radiation, so the spectrum will be different from that of the electromagnetic 
bremsstrahlung . 



II. FIELD EQUATIONS IN QUASILINEAR FORM 

Consider a system of point particles, interacting by non-gravitational fields (scalar ijj or massless vector A^) 
and moving in a self-consistent gravitational field described by the metric g^ v . The action can be presented as 
S = S p + jS^ + Sa + S g , where S p is the sum of particle actions including non-gravitational interaction terms 



S p = - J(m + fyj + eA^ds, 



SUand Sa are scalar and Maxwell field actions 



(1) 



(2) 



1 This is English translation of the Moscow State University preprint issued in Russian in 1980 The authors are grateful to 

Pavel Spirin for providing an electronic version for Fig. 1. This submission was supported by the RFBR grant 08-02-01398-a. 
* Electronic address: galtsov@physics.msu.ru 
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and the gravitational lagrangian is taken in the two-gamma form: 

s g = J cV^d% c = -^ J<r (r^r£ Q - r^r^) , K 2 = SnG. (3) 

Assuming gravitational field to be negligible at spatial infinity, we choose asymptotically Minkowskian metric 
r/ M „ in this region and introduce the (non-tensor) metric deviation variable 

h^v = g^u - Viiv ■ (4) 

Then for r — > oo one has /i M „ — > 0, but h^ v are not necessarily small everywhere. By convention, the indices of 
the quantities h^ v , and r/^ will be raised and lowered by the Minkowski metric r/^ , while the indices of 
the true tensors are operated with the metric g^. 
Introducing an antisymmetric tensor density 

jja.pX = g ^ g aX g 0u _ gaPgXv) _ (5) 

one can present Einstein equations in a divergence form 

(HJ 0X 9X» /y/=g) a = -2k 2 V~9 (TZ + t») , (6) 
where is the total matter stress-tensor, and t " is Einstein's canonical pseudotensor 

C - 4= 9f d{V ^ B £) - 6»C. (7) 
Maxwell equations can be written in a similar form 

(if J* 3 * g x , /y/=g) a = -^V~9f, (8) 

where the vector-current is 

f = Y j ejdsx» 5\x - x{ S ))/^—g. (9) 
Finally, the scalar wave equation reads 

{Anar V=9) iV = *kV=9t, (10) 

with the scalar current 

T = J2fJ ds 5 ^ x - < s ))/^9- (U) 
Particle equations of motion generically read 

^- [(to - ftp)} = ^g^^^x' 3 - mftp^ + eF^ x v . (12) 

All the above equations are exact and can be regarded as a system defining the particle motion and evolution 
of the scalar, vector and gravity fields and in a self-consistent way. However, since the notion of delta-functions 
is not defined in the full non-linear general relativity, we can deal with point particles only perturbatively, 
expanding all the quantities in formal series in the gravitational coupling k. For this one has to pass first to 
quasilinear form of the field equations. For Einstein equations we single out the linear part of the ff-tensor: 

jjav(3X _ ^afi^Xv _ ^aX^fiv _|_ £av@X _|_ j^avfiX ^3) 

where C au ^ x joins terms linear in h^ v : 

c *vf3X = ^olX^v + ^iv^X _ ^fi^Xv _ ^Xv^fi^ = ^ _ l^"^ (14) 

while U avfsx denotes non-linear in terms. Gravitational equations can be now presented as 

£ af) VftXVvT = 2k Tpi,, t^ v = T^u + S^v, (15) 

where in S^ v all the non-linear terms are collected. To calculate gravitational radiation one needs only terms 
mia.Hra.tir in h - 
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Maxwell equations can be rewritten similarly: 

A a , 0u (v afl V pl/ ~ V av V^) = Mf + Sn y/=g, (16) 
where an effective "gravitational" vector current is given by 

^—gS» = [(l/y/=g - 1) fa a V ~ »7 a V) + (C" af> + A^ Q/3 ) Ap, a ] v . (17) 
In the scalar case we obtain similarly: 

tZ = -4*f(T + S) , (18) 

where 

S = (i/4tt/>£ , <r" = (V^gg^ - vntv ■ 

This looks as the flat space wave equation for the spin zero field, with an important difference, however, that 
the "source" depends explicitly on tp. 

Now we can further simplify the quasilinear equations (which are still exact in all orders in k) imposing the 
gauge conditions 

^ v ,\r, uX = 0, A^iT = 0, (19) 
which are consistent with the field equations by virtue of the identities: 

Wf A =0, (V=5 5") ip = (V=5j'V (2°) 
In this gauge Einstein and Maxwell equations read 

□VV = 2k 2 t m „, nA„ = -<iny/=g~{j x + S X )t] X ^ (21) 

with □ = -dxd T 7] XT . 

III. SCALAR BREMSSTRAHLUNG UNDER GRAVITATIONAL SCATTERING 

Consider two point masses m\ and m.2, one of which (mi) carries the scalar charge /. Particles interact via 
gravity and the systems emits both gravitational and scalar waves. In this section we calculate scalar radiation. 
The action reads 

S = - J (mi + .M(x))^7^ds -m 2 J ^ 2 ds + V^9~d D xd^ (22) 

where x 2 = g a fjx a x^ , dot denotes differentiation with respect to the interval s, and the metric signature is 
mostly minus. 

The full system of equations describing the collision consists of Einstein equations, scalar field equation and 
particle equations as given in the previous section. The total loss of the four-momentum during the collision 
can be presented as 

ap£ = / at It 1 ; 1 do, , (23) 

J — oo J 

where 

is the energy- momentum tensor of the massless field tp, and integration is performed over the sphere of infinite 
radius. The non-zero contribution comes from the terms in , which fall off at infinity as r~ 2 . Also, without 
changing the integral, one can add to the integrand the total derivative over time (ip'^ip' — (l/2)7y A ' O '0 iCt V' Q ),o- 
Then, applying the Gauss theorem we can transform the Eq. (1231) to the following form: 



AP S =—\d i x W£ ■ (25) 



1 

To exclude an infinite self-energy part it is enough to substitute as t\i' v v the right hand side of the Eq. ([18 
while as ip - the t-odd part of the retarded potential 



M x ) = -1^X2 I d 4 ^-^e(fc°)^(fc 2 )[r(fc) + S(k)] , (26) 
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where e(k°) — 9(k°) — 9(—k°) is the sign function and the Fourier transforms are defined as 



T(k) 



d 4 xe lkx T(x) 



and similarly for S(k). We obtain: 

AP s = J~2 [ d 4 kk»6{k°)8{k 2 ) | r(fc) + S(k) \ 2 

Z7T J 



(27) 



(28) 



This expression is analogous to the usual one in electrodynamics, differing from it by presence of the non-local 
current S(k) which we will call the stress current. 

To find r(k) and S(k) we will solve the Einstein equation, the particle equations and the scalar field equation 
(IT51) expanding g^ v , ip and x^ in powers of the gravitational constant G. The actual expansion parameter in the 
ultrarelativistic collision problem will be the ratio of the gravitational radius of one of the particle to the impact 
parameter. One can show that approximation is valid if the particle scattering angle is small with respect to 
the radiation angle [13] 



G(mi + 1112)/ pv 2 < 1/7 , 

where 7 = (1 — v 2 )" 1 / 2 , v - the relative velocity of the colliding particles, p - the impact parameter. 
We parameterize the particles world lines as 

x^^ + ipl/m^sx + xtisx) , 

^2 = (P2/ TO 2)S2 +X 2 l {s 2 ) , 



(29) 



(30) 



with V— A 2 = p. (Api) = (Ap 2 ) = 0, A M - is the four-vector which in the rest frame of the second particle 
takes the form (0, p). Here and below we use brackets (. . .) to denote scalar products with respect to Minkowski 
metric. We choose the initial conditions 



z£(-oo) = dx£/ds(-oo) = 0, a = 1, 2 



(31) 



so that p£ is the four- momentum of the particle a before the collision, and rf 1 'PaPa — m a'i * s the correction 
due to the gravitational interaction. 

In the lowest order in gravitational interaction the correction to the space-time metric rj^ u due to the second 
particle reads 



K v {x) 



2G 

~^2 



u v 
P2P2 



1 



k l 



Substituting this into the equation of motion of the particle m\ we find 



~x»M 



G 



d i q 



Hqp2) 

q 2 {qpi) 2 



2(«Pi 



(PlP2)p- 



(PlP2) 



To calculate the Fourier-transform of S we use the expression for ip in the lowest order (zero order in G) 



tp(x) 



27T 2 



(32) 



(33) 



(34) 



The Fourier-transforms r(fc) and S(k) can be computed as follows. Using the integrals 



d4qe -. qA S(<lP2)5( q pi-k Pl ) = 

q l I 



d i qe 



-iqA 



6{qp 2 )S(qp 1 - kpi) ml{kpx) 

lf = Z7T 



/ 3 



fl (P1P2) IX 

Pi — P2 

mi 



K ( Zl ) -i(k Pl )A 



Zl 



(35) 



we obtain 



-r(fc) = -4Gmim2e Ife 



1 - 

z(fcA) 



(mim 2 ) 2 \ (P1P2) z 2 (mim 2 ) 3 



2/ 2 

•m\m 2 



I z x 2P 
{m 1 m 2 ) 2 \ Ki(zi) 



K (zi)- 



2P 



Zl 



(36) 



where I — {(P1P2) 2 ~ (mifr^) 2 } 1 / 2 , Z\ = \[— A 2 (fcpi)m 2 7 1 , z 2 = \/—A 2 (kp2)miI 1 Ko,i are the 

MacHnnalrl functions Lnrpnty.-i invariant intptrrals (IS fill an ho rnnvcnipntlv rnmniitprl in the rest, framp nf thp 
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second particle p 2 — (m2,0). Then two integrations from four are performed using the delta- functions, while 
the remaining two-dimensional integral in the plane orthogonal to pi is computed using the polar coordinates. 
The angular integral is standard, and the last one is done by contour integration. 

Integration over d 4 q in the expression for S(k) can be done using the Feynman parameterization. We obtain 

S(k) = 4GIe* kA z 2 2 I' dxe~ lx ^ , (37) 



o 



z(x) 



where 



z(x) = y/£ 2 (x); ^(x) = (1 - x)z x p%/m 2 + xz 2 p t ^/m 1 , 



z(0) =z x = V^£Z(k Pl )m 2 /I; z(l) = z 2 = V^(kp 2 )m 1 /I . 

The Eqs. (j3"6")l and (j3"T)) are obtained under the only restriction (f2l?|). they are valid for arbitrary velocities v. In 
the rest frame of the second particle p 2 = (m 2 , 0) we further specify the coordinates so that 

p = (0,0,p), p 1 = (0, Pl ,0), 

k = uj(sm 9 sin </>, cos 9, sin 9 cos 4>) . 

Consider the case of non-relativistic velocities. For «« 1, z(x) = z\ = z 2 = ujp/v, then the integral in (|37[) is 
easily done and we obtain 

= 2Gmim 2 ^ p sin g cos ^ _ . ^ & ^ 2 

= 4Gmim 2U aii: 1 (a)e- M ' ,siliecos ^ , (39) 

where a = cup/v. From (|38[) and (1391) it is clear that S'/r ~ w 3 , so in the lowest in v approximation radiation is 
entirely determined by the local current (f38|) . Substituting it into the Eq. (f28|). after some simple transformations 
we find 

d^Es = { Gfm 1 m 2)2a2[cos2 2 + ^ ^ 2 _ 
aw ail npv 

One can see that for small velocities the characteristic radiation frequency uj ~ u/p is inverse to the effective 
time of collision p/v. 

The total energy loss during the collision can be obtained integrating ([40]) over angles and the frequency: 

AE S = {ir/6)(fG mi m 2 ) 2 /p 3 v . (41) 

In the ultrarelativistic case (7 ^ 1) the effective spread of the stress current S(x) is of the order of p. So it can 
be expected that for the wavelengthes A ^ p the source with act as point-like. Indeed, for oj ^ p^ 1 (A > p) 
the argument of the Macdonald functions (|3"o| and (|37l) is small for all values of parameters, and with account 
for the leading terms we obtain 

AGm 1 m 2 sin cos S(k) . . 

r(k) = -i i-i *1 _LJ„ wp< i 42 

7 uipo 1 T(k) 

where 8 = 1 — v cos 9. Substituting (|42)) into ([28]) we find: 

d£ s 16 {fGm im2 f 2 _ a 

-7— = o 5 7 , w < p . (43) 

du) 3tt p z 

In the frequency region lo ^ p~ l the contributions of the local and the non-local currents are of the same order. 
In this case for the spectral distribution of the total emitted energy we obtain: 

w / / d Z dr l ,i I ^v,/,,, I _ 9 m/ 9 In 3 ■ ( 44 ) 



tt Jo Jo (1 + ^ 2 ) 3 / 2 (l + 77 2 ) 1 / 2 r? 2 

For p _1 < a; <C 7/p 

d£s . (/Gtoito 2 ) 2 2 fujp\ 2 

^7 = As ^ T — > (45) 
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where 



-Vs = f~ / <k»f" 8 In 3 + Vl + x 2 ) 

07T 



1) 



For relatively high frequencies ui ^> 7/p the integral in (|44|) is formed at £, 77 <C 1. So approximately 

d£ s „. 4(/Gm im2 ) 2 _ f /^\ ln *fW e -^ (46) 



rfw p 2 V 7 / 7 

The expressions (|43p . (|45)) and (|46|) together describe the behavior of the spectral curve. It follows, in 
particular, that in the spectral distribution there is a maximum around the frequency uj ~ 7/p- The total 
energy loss during the collision is obtained integrating (|4"4"]) over the frequency: 



(/Gmim 2 ) 2 3 . 3G 77 
"7 7,A S = — + - 



AEs = As ^^l 7 3 , A S = ^ + ^ - 2vr ~ 1.51 , (47) 



where G is the Catalan constant. 

Let us compare these results with the case of the electromagnetic interaction in Minkowski space, when the 
source term in the equation for tp does not contain the stress current S(x). Suppose that the colliding particles 
are electrically charged (ei,e2) and neglect gravitational interaction with respect to electromagnetic one. Then 
as the source r(k) in ([25)) one has to use the Fourier-transform of the trace of the particles energy-momentum 
tensor. After similar calculations we obtain: 



T ^ _ 2e 1 e 2 (m 1 m 2 f j kA 



(P1P2) - m 



1 (fcpi) 



K ( Zl ) - i(kA)(p lP2 )^^ } . (4S) 



In the ultrarelativistic case (7 3> 1) the spectral- angular distribution is dominated by the second term in 
In the rest frame of the second particle we find in the leading order in 7: 

^.•^./•^(..iW, ,49, 
dui up 1 J z x \ x J 

where z = wp/27 2 . Using the asymptotic expansions for the Macdonald functions for small and large arguments, 
from we find for u <C 7 2 /p: 

dE s _ 4 (Me 2 ) 2 

1 — ~ o 2 ' ^ ' 



while for high frequencies co 3> 7 2 /p 

^ = 2 (/^)V e -?, (51) 

dui p z ujp 

Note, that for the local source case our methods gives the energy loss without restrictions on the relative 
velocity of collision. Indeed, substituting (|48p into (f^5|) and integrating over frequencies and angles we obtain 

7T (/e x e 2 ) 2 / 2 1 

— 3 — 7 + q 

tip 13 \ 3 



AEs = ^-^^(^ + ^) . (52) 



Thus we see that there is substantial difference between the spectrum of the bremsstrahlung from gravitational 
scattering and that in the case of electromagnetic interaction. In the first case there is a maximum at w ~ 7/p, 
while the spectral distribution (|49|) is monotonous function of the frequency. For gravitational interaction 
the exponential cut off corresponds to the frequency oj ^ 7/p, and not to u ^ 7 2 /p as in the case of the 
electromagnetic scattering. Finally, the total energy loss at gravitational scattering ([4"Tf is 7 times less that the 
corresponding quantity in the electromagnetic case (|52|) for the same scattering angle, i.e. under the condition 
G7712TO17 ~ eie 2 - 

These properties can be qualitatively explained by the presence of the non-local (in terms of the flat space-time 
picture) stress-current source in the equation for the radiated field ip in the case of gravitational scattering. This 
current has and effective transverse dimension of the order of p and longitudinal of the order of p/7 (7 times 
smaller due to the Lorentz contraction). For large wavelengthes (A 3> p) the source non- locality is insignificant 
and the low frequency limit is the same as for the electromagnetic interaction case, when there is no non-local 
term at all. For A < p radiation from the most distant elements of the source exhibit a destructive interference 
for the angles close to 7r/2, which leads to the gap in the spectrum. Finally, for A < p/7 the conditions for 
destructive interference are fulfilled for the forward direction, in which the most of the energy is emitted. This 

leads tn substantial decrease nf the radiation 
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IV. ELECTROMAGNETIC BREMSSTRAHLUNG UNDER GRAVITATIONAL SCATTERING 

The case of the electromagnetic interaction is rather similar. Let the particle mi carry the electric charge e\. 
Using analysis of the Sec. 2 we can present Maxwell equations as follows: 

{rirrf p F a p\ v = -4tt(J^ + S*) , (53) 



where the stress-current is 



S" = off , (54) 



/fir** 
d8=±5(x - Xi(8)) 
ds 



(55) 



Imposing the flat space Lorentz gauge on the four-potential A^: 

rTA M , v = , (56) 
we cast Maxwell equations into the form convenient for iterative solution: 

rTv^^fi = 4^(J" + S") . (57) 
It is convenient to choose two linearly independent polarization vectors as 

e» = X+eTP'KpipVu, , e£ = \ge^k v e^ p p 2a , (58) 

A = (-P 2 )-V2 ; p m = {kp2)p » _ {kpi)p » _ Xg = {kp2 yi . (59) 
They satisfy the following conditions: 

(ege^) = (keg) = [ke^j = , (e^e^) = (egeg) = -1 (60) 

and in the rest frame of the second particle read: eg = (0, eg), = (0, e^) , where eg and are unit vectors 
along 9 and (f>. 

The expression for the four-momentum loss due to electromagnetic interaction with polarization A (A = 9, <p) 
can be derived analogously to the Eq. (|2"5]l and reads: 

AP e^ = ^2 / d 4 kk^(k°)6(k 2 ) I I^(k) I 2 , (61) 

where I^(k) — ?7 Q ^e"(J' 3 (fc) + S l3 (k)). As in the scalar case, one has to retain in only terms falling off 
asymptotically as r~ 2 . In this approximation 

S»(x) = -(1/4tt)(F£/i ct1/ + F v a U ia - {h a j2)F^ v ), v . (62) 

The subsequent calculations are similar to the scalar case. The Fourier-transforms of the currents (|55|) and (|6"2"T) 
are computed in the full analogy with the previous section resulting in 

JM(fc) ^f^i^i 1 - ^^)*>(*0[(**0tf - OtoOrfR 

+ ^(l + ^f^)^(, 1 )PAK-(^ 1 )A^]} , (63) 

4Gm 2 ikA f 1 j -ix(kA){ A 2\( k P2) 



S»(k)=-^p-e ikA / dxe~™^ 







x rr- 1 



(1 - a;)( fcpi )^ lP2 ) 



2/ 2 



2/ 2 

Ki(z(x)) ( M A . (pipa) / , 



x ((fcp 2 )A^ - (fcA)p^ - ijf (*(a:)) f (*A)p? - (Api)A") } . (64) 



Tn (lfi4D terms, Tvronortinna.l to k} 1 are omitted since thev Ho not cnntririiite to radiation hv virtue of (IfiOL 



8 



For small relative velocity (v -C 1) radiation is generated predominantly by the local current J M , since 
S x /J x — v 2 . In this case 

d 2 E e em _ ( e i Gm 2Y shl 2 Q o 2 [x g( ) + ctgtesm* <PK 2 (a)} 



dujdfl \ npv 



dutdfl \ npv 
Integrating over frequencies and angles we get 

g _ 7ir {e 1 Gm 2 f 
em 48 vp 3 



- ' 1 cos 2 4>a 2 K( (a) . (65) 



AEL= ^(j^f. (66) 
lb vp 6 

For ultrarelativistic collisions (7 > 1) in the low frequency range {uj <C 1/p) contribution of the non-local 
stress current is relatively small, S x / ' J x ~ up <C 1, and we have: 

dE e em 8 ( ei Gm 27 ) 2 



dcu 37T p 2 
dEf m _ 8 ( ei Gm 27 ) 2 



g?u; 7T p 2 

For w ^ in the leading in 7 approximation the spectral distribution of the radiated energy is given by 



(67) 



«w I" Jo Jo 



i+2e \+e+v 2 

■ In 



2 



( 1+ £3)3/3( 1 + J7 2)l/2 ^ 

In (|68p we performed summation over polarizations. 

Using (|68|) one can show that for p^ 1 < 10 <C j/p the spectral distribution behaves as follows: 

d J^ = (e lGm2 ) 2 (l) 2 (^)l^, (69) 
du \pj \ 7 ) up 

while for the frequencies ui ^ 7/p 

^ = 4 (ei Gm 2 ) 2 f lV f ^ In ^l e -^ . (70 ) 
\pj V 7 / 7 

Comparing ([BT| . (f6"!)f and ([70)1 one can notice the fall off in the spectrum in the frequency range u ~ p^ 1 and 
the maximum at w ~ 7/p (Fig. 1). 
For the total energy loss we obtain: 

A£? em = A em - - — , (71) 

where A em = 5G/2 + 43/12 - 7r ps 2.75. Splitting on polarizations is given by A em A x m (A e em « 1.75, Af m rj 
1.00). Our result (|71l) qualitatively agrees with that of [li| but differs from that of [22j by absence of the factor 
In 27. 

In the case of both particles electrically charged with large charge to mass ratio in geometric units one can 
neglect gravitational interaction and the bremsstrahlung problem is simplified considerably. Then the stress- 
current = 0, and the radiation amplitude is fully given by the local current. Consider for simplicity the case 
rn 2 3> mi. Then the Fourier-transform of the current is given by 

J"(*) = 2(eie2)( r iW2)2 ^ A (^ - T^rfWo + i^((kpi)A» - (^AK)^M) . (72) 
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In the non-relativistic case (v <C 1) the spectral-angular distribution of the emitted energy, as it can be expected, 
is given by the Eq. (I65|) with the substitution Gm\mi — > &\&i- As before, two independent polarization states 
are given by the unit vectors (|58p. (I5T)1) . Using the Eqs. (|72p. ([5"8"]) and ([S"9"]) and passing to the rest frame of 
mi one finds with account for (1611) the following expression for the energy loss due to electromagnetic radiation 
with the polarization A (A = 9,(j)): 

AE 9 = Hi±ll( 1 2 + l 
em 64 mjpZv [J 3 

™L = + 1)- (73) 



64 m\p z v \ 3 

Note that Eqs. (|75|) are valid for an arbitrary relative velocity of the colliding particles. In the ultrarelativistic 
case (7> 1) the spectral distribution is given by the second term in (|72[) . so in the leading approximation in 7 

^_i^,r (b ( 1 _2» + 2£W (74) 

duo it m\p z J z \ x ar J 
where z = iop/2^ 2 . In the low-frequency limit 00 <C 7 2 // the Eq. [THhas the form 

dE em _ 8 (e 2 e 2 ) 2 



dui 37r m\ p 2 



(75) 



which coincides with (|67p. if both results are expressed in terms of the scattering angle. At high frequencies 
lo 3> 7 2 //3 the spectral distribution has exponential cut off: 

d^n = (flfafe-^ . (76) 
duo mfp z 

The numerical curve for the spectral distribution is given in Fig. 1. 

Comparing the Eqs. (J67J), (JSU), ([70]) and ([71]) with the Eqs. ([73]), ([75]) and ([76]) one can see that the difference 
between spectral properties of radiation for particles interacting by gravity and by non-gravitational forces is 
similar for scalar and electromagnetic radiation. 

V. GRAVITATIONAL BREMSSTRAHLUNG 

Consider now the system of two gravitating point particles m\ and mi. We choose coordinates in such a way 
that the metric perturbations be small at infinity when particles are at finite distance from each other. Then 
we can treat the particles at t = ±00 as free and the metric to be flat (excluding the self-field of each particle 
in its vicinity which can be removed by classical renormalization). Denote the covariant components of the 
4- moment a as 

p a = lim m a u a , u a = g^dx^/ds , p' a = lim m a u'°" , a = 1,2 . (77) 

The change of the total four-momentum of the system is due to radiation friction acting on the particles. 
Although for two relativistic gravitationally interacting particles it is problematic to find the gauge independent 
local reaction force, one can still find in a coordinate independent way the expression for the total momentum 
loss during the whole collision time: 

, f°° du a 



o=l,2 Q=-l,2 



This quantity can be shown to be independent on the coordinate choice if the coordinate transformation preserve 
the above asymptotic conditions. Using the equations of motion we find 



AP M = i / d^xV^ggv^T™ . (79) 



2 

Since the covariant derivative of the stress-tensor is zero T^.jy — 0, wc hcLVG 

gu*,^gT v ° = 2(V=gTp, u . (so) 

Now we make use of the conservation equation 

u/=flfT.?+tf:vi =0, (8D 
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where t" is the Einstein pseudotensor. In our approximation it will be enough to keep only quadratic terms in 
h„ 



1 



32ttG 



As a result, we transform the momentum loss to the form 



AP„ = - J d 4 x(V=9 ty, v . 



We assume the gauge ipff — and calculate the divergence of the pseudotensor to get 

1 



AP„ 



32ttG 



d A xh afi ^<f> v . 



(82) 



(83) 



One can show that the lowest order giving non-zero contribution is the second (or the first post-linear order). 
Using the Einstein equations in quasilinear form, as given in the second section, we perform transformations 
similarly to the electromagnetic case introducing the polarization tensors for gravitational waves. The final 
expression for the loss of the four-momentum on gravitational radiation with the polarization A reads: 



G 



APW" = / d A kW6(k )5(k 2 ) | T< A >(fc) 



(84) 



where = e$T» u , r^ u = + , 

-16wGS^=h^h u ^ a - h a /h^ - {l/2)h^h a p, u + 

- (l/2)Ji£V + (V2)^( 2 ^^,A " Kp,xh aX ^ + {3/2)h a p, x h a ^ x ) , 



(85) 



and it is assumed that all contractions over indices in (|84f and (|85f are performed with Minkowski metric ij^ 
It is convenient to choose as two independent polarization vectors the quantities 



e* = (1/V2)(eiet + e'e*), e+ = (l/V2)(e!eS - e*e% 



(86) 



-~fJ,U Z 



- 0) e nu — e vui ^ L e x — 0, A — X, 



(87) 



The subsequent calculations are essentially similar (though more lengthy) as for the scalar and electromagnetic 
radiation, so we give the final result. The amplitudes T^ v and S^ v in an arbitrary Lorentz frame read 



T^{k)=T? v {k)+T% v (k) , 



where 



Tr(fc) =AGe 



- 1 



_(■,_ {mim 2 ) 2 \ (gipg) m 2 z 2 _ { 
V 2P ) I mi zi 

(m 1 m 2 ) 2 \ (pipa) 



(raim2) 2 \ 77^2 



2/ 2 



2I 2 



(m 1 m 2 ) 2 \ml K\{z\) „ „ 
-PiPr 



»(*pi) 1 + 



(mim 2 ) 2 \ m 2 K\(z\) 



2/ 2 



I 



2/2 



Zl 



Tj™ (lf)2), 



S^{k) = AGIe lkA J dxe-* x(kA) l(l + 
o ^ 



z{x)K 1 {z{x)) 



2/2 



1)17 



m 1 m 2 {p 1 p 2 ) ( (p\Pi) . f1 s 

1 z 2 x + zi(l - x) 

l z \mim 2 



(mim 2 ) 4 / (pip 2 ) 

2/ 4 \ TYl\TYl 2 



(m 1 m 2 ) 2 \m 2 2 „ „ 

It7^o(zWKpi+ 



2/2 



z 2 x + z\(l — x) 



Ki(z(x)) 

PiPi 



z(x) 



- IZ 2 - 



m 2 



m\m 2 \ ll A J \mim 2 J 



K {z(x))[p>i 



A" 



Pv 



A^ 



(89) 
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Note that in the electromagnetic case the local and non-local currents are separately gauge invariant, while in 
the gravitational case only the sum T^ v + T£ v + is independent on the gauge choice. This allows to change 
contribution from separate terms choosing suitable gauge. In particular, in the gauge (|86p the contribution from 
T^ v is zero. The subsequent calculation will be performed in the rest frame of the second mass p 2 = (m 2 , 0). 

In the non-relativistic limit (v <C 1) z(x) ~ z\ w z 2 w ujp/v , so the integral (|89p can be easily computed. 
The contributions from (|88[l and f|89[) turn out to be of the same order, and taking into account (|84|) . one finds: 

d 2 E+ = G 3( 2 )2 ^ 2 ^ ^ ^ ^ + Ko{a)] 2 + 

audi I 7r z /cr 

+ [sin 2 6K Q {a) + (sin 2 6 + sin 2 <j) - cos 2 cos 2 0)aJfi (a)] 2 } , (90) 



cP^fc = 4 G3(TO 9 1 ^ 2)2 a 2 {cos 2 gcos 2 c/m 2 K 2 (a) + sin 2 9 sin 2 0[Jfi(o) + ftA' (a)] 2 } , (91) 
dujdll TT z p z 

and after the integration 

+ 43 2 77T G 3 ( mi m 2 ) 2 v 343tt 

The Eqs. (jnnj) and (jH]) coincide with those given in fl7j |. 

An ultrarelativistic case is considered similarly to the previous sections. For u < wc obtain an expression 
coinciding with the result of application of the low frequency theorems [27] 



dE+ 32G 3 ( mim 2 ) 2 f 7 



duo 'in \p 



(93) 



^_64 G 3 (mim2)V7V ^ n27 _r 



Note different dependence of (|93| and (|94|) on the energy. 

For the frequencies uj ^ p^ 1 the leading in 7 approximation the spectral distribution of the gravitational 
bremsstrahlung summed up over polarizations is given by 

00 00 

dE,, = iggw^^ / / ^-^yj^^P 1 + 8 f + ^ 1/9 m . (95) 



7T J J (l+e 2 ) 3/2 (l + ?? 2 ) 1/2 ^ 



In contrast to the previous cases, the spectral curve is monotonous function of the frequency, and for relatively 
small frequencies p^ 1 < lj -C 7/ p the spectrum falls off logarithmically 

dE gr _ 64G 3 (mim 2 ) 2 M 2 2 7 



duj 7r V /° 7 eP up 

while for w 3> 7/p - the fall off is exponential 

dE gr _ 3 2 fj\ 2 fujp\ 4e c a;p^2^ £ 



- InV- ' ( 96 ) 



4G a (m 1 m 2 ) 2 -U ^ In ^^e"^ . (97) 
VP/ V 7 / 7 

Note that for cj = p -1 (|96| with logarithmic accuracy coincides with (|94|). while (|96|) and (|97|) by the order 
of magnitude are compatible at cwy/p. So the Eqs. (|94|) . (|96|) and (|97]l together covers the whole frequency 
spectrum. The total radiated energy is obtained integrating (|9"5j) over frequencies 

AE gr = A 9r G 3 (TOim 2 ) 2 ( 7 /p) 3 , 



A gr = vr + 35G/2 -211/12 29 . (98) 

The result (|9"5|) by the order of magnitude coincides the the results of refs. HI]) but differ from [22], [25[ 

by the absence of the factor ln(27). The frequency distributions of scalar, electromagnetic and gravitational 
radiation under ultrarelativistic gravitational scattering are shown in Fig. 1. 

Consider now for comparison gravitational radiation under collision mediated by non-gravitational forces. Let 
both particles be charged (with e\ and e 2 correspondingly) with large charge to mass ratio, so their gravitational 

interaction can be neclecterL Then in the lowest in d a.nnroximation n.... = n = ft and r.... = T....4-&T..... 
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dm 




Figure 1: The spectral distribution of scalar (green), electromagnetic (red) and gravitational (black) radiation under 
gravitational scattering for 7 = 1000. 



where fT^ is the energy- momentum tensor of the electromagnetic field. Calculations shows that with the same 
accuracy 



T[ lu {k) 



{pip 2 ) + m 1 m 2 — 1^0(2:1)^1^1 +mfi^ (2i)(PiP2 + P1P2)- 
z\ J 



Zl 



Zl 



T£"(k) = e lfcA Tf^(l o 2) , 



? T^(k) = 2eie 2 e lfeA / dxe'^^ 



1 / O1P2) 



z(x)Ki(z(x)) 



A M A 1- 



(mim 2 ) 2 xxPiK 
-A (z(x)) — — 



/ 2 



(mim 2 ) 2 / (pip 2 ) 



/ 2 
,mim 2 



TOim 2 
O1P2) / (P1P2) 



-z 2 x + Zl(l — x) 



ni 2 



z 2 x + 2i(l — z) 



(P1P2) / (P1P2) 
J \mim 2 
I 



z 2 x + zi(l — x) 



m\m 2 



n 
■22 



A'i tfpg 



z(x) m 2 



2m\m 2 



z 2 



K 



A^ p\ A v p^ 



V-A 2 mi V-A 2 mi 



(100) 



and in the chosen gauge the contribution form T^ v is zero. For 7 3> 1 in the frequency region w < one 
obtains the results (|96l97p , in which the gravitational deflection angle should be replaced by the electromagnetic 
one 4Gmim 2 7 — > 2eie 2 . 

For uj p^ 1 the non-locality of the source due to presence of the term pT^ leads to destructive interference. 
Like in the above cases of the gravitational interaction the non-local source f T^ v gives rise to two impulses of 
different duration, one of which comes to the observation point in the counterphase with the one due to T^ v ' . 
But contrary to the case of gravitational interaction, the contribution of T^ u is canceled only partially. As a 
result, the spectral-angular distribution for uj ^ p _1 and 7 lwill be 



d 2 E n 



2G(eie 2 ) 



duidcosO 

In the leading in 7 order we obtain 



uj 2 Kf(ujpd) sin 2 6 [ 1 



sin 2 9 
2^fS 2 



dE gr 4G(eie 2 ) 5 



z / dxK l (x) 



Az 2z 2 



dui up 

which has the following asymptotic behavior for p _1 < uj <C 7 2 /p 

dE, 



-3 



duj 



v _ 4G(eie 2 ) 2 ^ 4 7 



e c uip 



for uj 3> j 2 /p 



dE, 

rl,., 



G(eie 2 ) 2 



1 2 \ -H 
— I e < 



(101) 



(102) 



(103) 



(104) 
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The total energy loss is 

AE gr = (7r/4) 7 2 G(e ie2 ) 2 /p 3 ■ (105) 

The result (I105P coincides with that of [2lj |. 

The spectral distribution of gravitational radiation in two considered cases has the following distinctive 
features. For it weakly depends on frequency and for the fixed deviation angle depends on the energy 

as 7 2 ln2 7 . For the frequencies p^ 1 < lu <C (jJ C t the spectrum falls off logarithmically, while for lj 3> w C r 
exponentially. But if for the electromagnetic interaction uj^ n — 2-f 2 /p, for gravitational interaction uj 9 J t = j/p. 
Also, for the same scattering angle the total radiated energy for electromagnetic interaction the radiative loss 
is 7 larger that for gravitational interaction. 



VI. LOW FREQUENCY LIMIT 



For u) — > the spectral distribution does not depend on frequency and one could hope to get a correct estimate 
for the energy loss under collision multiplying the Eq. (|67|) or (|93|) and (|94|) on a suitable frequency cutoff. For 
radiation of the point particle in the flat space (in the case of non-gravitational interaction) the cutoff frequency 
in the classical spectrum is estimated kinematically as an inverse time of the formation of radiation in the given 
direction and it is given by w^™ ~ j 2 /p. In the gravitational case a similar estimate is wf, r ~ 7/p which is 
confirmed by an accurate calculation. Now, it can be easily seen that the low frequency approximation gives a 
correct estimate of the total radiated energy in the electromagnetic case, but gives an wrong factor In 2 7 in the 
gravitational case. The reason of this discrepancy lies in the fact that the fall-off in the spectral distribution 
in the gravitational case corresponds not to the frequency uj cri as it is assumed in the low-frequency approach, 
but to w ~ p _1 (see (I96p and (|103p V Logarithmic fall-off in the high frequency region u> ^ p" 1 cancels an extra 
logarithmic factor. This explains the difference between our result (|98p and that of [24l l25j|. 



VII. METHOD OF VIRTUAL GRAVITONS 



The spectral density of the wave packet of equivalent gravitons imitating the gravitational field of the ultra- 
rclativistic particle is given by 



I gr (u,p) = G{m 2 /irp) 2 (ujp/j) 2 Ki(ujp/'y) 



(106) 



(this results differs from that of the ref. [22j ] by a numerical factor). The spectrum diverges for ui — > 0. 
This means that it can be applied only for sufficiently high frequencies. Applying this spectrum to compute 
bremsstrahlung (electromagnetic or gravitational) under scattering of the fast particle on the fixed center one 
has to introduce the frequency cutoff. The results differ from those obtained in this paper by a factor In 27. 
Thus, contrary to the electromagnetic case, where method of virtual quanta gives the correct answer in the 
ultrarelativistic limit, in the gravitational case this method fails. The reason is that the spectrum of virtual 
gravitons describes correctly the frequency range lj 3> 7/p, which, as we have seen, is negligible in the total 
radiation due to non-locality of the effective radiation sources. Indeed, let us consider radiation in the forward 
direction. For 8 = the integral over the Feynman parameter can be computed exactly and we obtain 



d 2 E a 



duidfl 



G 3 (mirri2) 

7T 2 7 2 



Ki(z 2 ) 



(107) 



At the same time the equivalent gravitons approach gives 



d 2 E, 



duodVl 



G 3 (mirri2) 

7T 2 7 2 



-to 2 K 2 ( Zl ) 



(108) 



One can see that the expressions (|107[) and (|108[) are compatible only for ui ^> -f/p. 



VIII. CONCLUSIONS 



We have presented Lorentz-covariant perturbation approach in General Relativity using the momentum space 
formulation similar to quantum field theory perturbation theory. The method consists in solving particles equa- 
tions of motion and the field equations iteratively in terms of the gravitational coupling constant. Gravitational 
radiation arises in the second order approximation. In terms of the flat space metric the source of the D'Alcmbcrt 
equation for the second order metric perturbation is non-local and contains the contribution from gravitational 
stresses computed in the first order. This non-locality results in 7 times lower frequency cutoff as compared 

tn thp rasp of nnn-ora.vitatinnal intm-artinn Fnr this rnasnn flip mpthnH of virtual pravitnns is not annlicahlp 
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for gravitational scattering of ultrarelativistic particles. The total energy loss in the rest frame of one of the 
particles is proportional to the third order of 7. Radiation from two colliding bodies looks as a collective effect, 
contributions form each of them can not be separated in a gauge independent way. 
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